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Abstract 

A study of Jjili-K^ and — > '(/;(25)7r* decays is performed with data corresponding to 

0.37 fb~^ of proton-proton collisions at = 7TeV. Their branching fractions are found to be 

J/i/'TT*) = (3.88 ±0T1± 0.15) X 10"^ and 
BiyB'^ ^ i){2S)-K^) = (2.52 ± 0.26 ± 0.15) X 10"^ 

where the first uncertainty is related to the statistical size of the sample and the second quantifies 
systematic effects. The measured CP asymmetries in these modes are 

A'I^p^ = 0.005 ±0.027 ±0.011 and 

^^(2S)^ = 0.048 ± 0.090 ± 0.011 

with no evidence of direct CP violation seen. 
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The Cabibbo-supprcsscd decay — > ^/;7r+, where 
represents either a J/ip or tp(2S), proceeds via a 
b — ^ ccd quark transition. Its branching fraction is 
expected to be about 5% of the favoured b ccs 
mode, B'^ — > tpK^ (charge conjugation is iniphed 
unless otherwise stated). The Standard Model 
predicts that for h — >■ ccs decays the tree and 
penguin contributions have the same weak phase 
and thus no direct CP violation is expected in 
5+ —5- ipK'^ . For B^ ipTT'^, the tree and pen- 
guin contributions have different phases and CP 
asymmetries at the per mille level may occur [I]. 
An additional asymmetry may be generated, at 
the percent level, from long-distance rescattering, 
particularly from decays that have the same quark 
content {D^D-, D*- , ...) g]. Any asymmetry 
larger than this would be of significant interest. 



In this paper, the CP asymmetries 



B{B- iIjtt-) ~ B{B+ V'7r+) 



B{B- i'Ti-) + B{B+ tjm+) 
and charge-averaged ratios of branching fractions 



1) 



B{B^ -^i'TT^) 

B{B^ iIjK^) 



are measured with the V reconstructed in the 
fX^ IJL~ final state. From the latter, B{B^ — >■ ^-k^) 
may be deduced using the established i?^ — > ipK"^ 
branching fractions [3]. The CP asymmetry for 
B+ ij{2S)K+ is also reported. B+ J/iPK+ 
acts as a control mode in the asymmetry analysis 
because it is well measured and no CP violation 
is observed [3]. Previous measurements of the 
5+ — >■ J/'07r+ branching fractions and CP asym- 
metries [HIS] have an accuracy of about 10%. The 
— > -0(2S')/i+ {h = K,tt) system is less precisely 
known due to a factor ten lower branching fraction 
to the hfifi final state. The world average for 
^V(2S)K -0.025 ± 0.024 [3| and there has been 
one measurement of ^'^(^■S)'^ = 0.022 ± 0.086 Hj. 

The LHCb experiment [7] takes advantage of 
the high bb and cc cross sections at the Large 
Hadron Collider to record unprecedented samples 
of heavy hadron decays. It instruments the pseu- 
dorapidity range 2 < 77 < 5 of the proton-proton 
(pp) collisions with a dipole magnet and a tracking 
system which achieves a momentum resolution of 
0.4 - 0.6% in the range 5 - 100 GeV/c. The dipole 
magnet can be operated in either polarity and this 
feature is used to reduce systematic effects due 
to detector asymmetries. In the sample analysed 



here, 55% of data was taken with one polarity, 
45% with the other. 

The pp collisions take place inside a silicon-strip 
vertex detector which has active material 8 mm 
from the beam line. It provides measurements 
of track impact parameters (IP) with respect to 
primary collision vertices (PV) and precise recon- 
struction of secondary i?+ vertices. Downstream 
muon stations identify muons by their penetration 
through layers of iron shielding. Charged particle 
identification (PID) is realised using ring-imaging 
Cherenkov (RICH) detectors with three radiators: 
aerogel, C4F10 and CF4. Events with a high 
transverse energy cluster in calorimeters or a 
high transverse momentum (pt) muon activate 
a hardware trigger. About 1 MHz of such events 
are passed to a software-implemented high level 
trigger, which retains about 3 kHz. 

The analysis is performed using 0.37 fb~^ 
of data recorded by LHCb in the first half of 
2011. The decay chain B+ iph+ , tfj 
is reconstructed from good quality tracks which 
have a track- fit P^r degree of freedom < 5. 
The muons are required to have momentum, 
p > 3GeV/c, and pt > 0.5GeV/c. Selected 
hadrons have p > 5 GeV/ c and > 1 GeV/c. 
The two muon candidates are used to form a ?/> 
resonance with vertex-fit < 10- The dimuon 
invariant mass is required to be within ^49 MeV/c^ 
of the nominal ip mass [3]; the asymmetric limits 
allow for a radiative tail. 

The reconstructed B^ candidate vertex is 
required to be of good quality with a vertex-fit 

< 10. It is ensured to originate from a PV by 
requiring Xip < 25 where the x^ considers the 
imccrtainty on track IP and the PV position. In 
addition, the angle between the B~^ momentum 
vector and its direction of flight from the PV 
must be < 32 (10) mrad for ^{2S)h+ iJ/iph+). 
Furthermore, neither the muons nor the hadron 
track may point back to any primary vertex with 
Xip < 4. It is required that the hardware trigger 
accepted a muon from the B~^ candidate or by 
activity in the rest of the event. Hardware-trigger 
decisions based on the hadron are neglected to 
remove dependence on the correct emulation of 
the calorimeter's response to pions and kaons. 

The candidates are refitted [8] requiring all 
three tracks to originate from the same point in 



1 



space and the ip candidates to have their nominal 
mass [3] . Candidates for which one muon gives rise 
to two tracks in the reconstruction, one of which is 
then assumed to be the hadron, form an artificial 
peaking background in the '0(2S')/i+ analysis. 
These candidates peak in the invariant mass distri- 
bution of the same-sign muon-pion combination at 
TO/^TT ^ 245 MeV/c^, i.e. the sum of the muon and 
pion rest masses. Requiring m^jr > 300 MeV/c^ 
removes this background. In 2% of events two 
candidates are found. If they decay within 2 
mm of each other the candidate with the poorest 
quality vertex is removed; otherwise both are kept. 

When selecting J/V'/i'^ candidates, a require- 
ment is made on the decay angle of the charged 
hadron as measured in the rest frame of the i?+ 
with respect to the 5+ trajectory in the laboratory 
frame, cos(0^) < 0. This requires the hadron to 
have flown counter to the trajectory of the i?+ 
candidate, hence lowering its average momentum 
in the laboratory frame. At lower momentum, 
the pion-kaon mass difference provides sufficient 
separation in the i?+ invariant mass distribution, 
as shown in Fig.[T] In the — > ■0(2S')/i+ analysis, 
the average momentum of the hadrons is lower, so 
such a cut is unecessary to separate the two modes. 




5000 5100 5200 5300 5400 5500 

m(J/\\m^) (MeV/c^) 

Figure 1: Distribution of cos(6'j^) versus the invari- 
ant mass of J/-07r+ candidates. The curved 
structure contains misidentified — !> J/if^K^ de- 
cays which separate from the — > J/-0 7r+ vertical 
band for cos (6*,*) < 0. The partially reconstructed 
background, B J/-4>Kn enters top left. 



Particle identification information is quan- 
tified as differences between the logarithm of 
likelihoods, IwCh, under five mass hypotheses, 
h G {tt, K, p, e, fi}. Separation of ip?:'^ candi- 



dates from tpK^ is ensured by requiring that the 
hadron track satisfies In Ck — In/^Tr = DLLktt < 6. 
This value is chosen to ensure that most 95%) 
— > -07r+ decays are reconstructed as such. 
These events form the "pion-like" sample, as op- 
posed to the kaon-like events satisfying DLLa'tt > 6 
that are reconstructed under the ipK^ hypothesis. 

The selected data are partitioned by mag- 
net polarity, charge and DIA^ktt of the hadron 
track. By keeping the two magnet polarity 
samples separate, residual detection asymmetries 
between the left and right sides of the detector 
can be evaluated and hence factor out. Event 
yields are extracted by performing an unbinned, 
maximum-likelihood fit simultaneously to the eight 
distributions of B invariant mass in the range 
5000 < TUB < 5780 MeV/c^ ^. Figure [H shows 
this fit to the data for i?+ — > J/tph'^, summed 
over magnet polarity. The 5+ — ■!/'(2S')/i+ data is 
shown in Fig. [3] 

The probability density function (PDF) used 
to describe these distributions has several compo- 
nents. The correctly reconstructed, B~^ — >■ tjjh'^ 
events are modelled by the function, 

which describes an asymmetric peak of mean fi and 
width cr, and where aL{x < fi) and aj^^x > fi) pa- 
rameterise the tails. The mean is required to be the 
same for ipK^ and tpir'^ though it can vary across 
the four charge x polarity subsamples to account for 
different misalignment effects. Table [T] shows the 
fitted values of the common tail parameters and 
the widths of the _B+ ~> ijjh^ peaks averaged over 
the subsamples. 

Table 1: Signal shape parameters from the B^ 
iPh^ fits. 





J/4, 


4,{2S) 




(MeV/c^) 
(MeV/c2) 


7.84 ±0.04 
8.58 ±0.27 


6.02 ±0.08 
6.12 ±0.75 






0.12 ±0.03 


0.14 ±0.01 


aR 




0.10 ±0.03 


0.13 ±0.01 



The misidentified •ipK'^ events form a displaced 
peaking structure to the left of the -(/'7r+ signal 
and tapers to lower mass. This is modelled by a 
Crystal Ball function [TU] which is found to be a 
suitable effective PDF. Its yield is added to that 
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Figure 2: Distributions of J/^h^ invariant mass, overlain by the total fitted PDF (thin line). 

Pion-like events, with DLLif^r < 6 arc reconstructed as J/i/jn^ and enter in the top plots. All other 
events arc reconstructed as J/ipK^ and are shown in the bottom plots on a logarithmic scale. B~ decays 
are shown on the left, i?+ on the right. The dark [red] curve shows the J/ifin^ component, the 

light [green] curve represents B^ — > J/^K^. The partially reconstructed contributions are shaded. In 
the lower plots these are visualised with a dark (light) shade for (i?+ or i?°) decays. In the top plots 
the shaded component are contributions from B J/ipK^n (dark) and B J/^pn^TT (light). 
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Figure 3: Distributions of B^ — > ip(2S)h^ invariant mass. See the caption of Fig. [2] for details. The 
partially reconstructed background in the pion-like sample is present but negligible yields are found. 
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of the correctly identified events to calculate the 
total number of ^K'^ events. 

The PDF modelling the small component of 
ipTT'^ decays with DLL/^^ > 6 is fixed entirely 
from simulation. It contributes negligibly to the 
total likelihood so the yield must be fixed with 
respect to that of correctly identified -i/jTr"*" events. 
The efficiency of the PID cut is estimated using 
samples of pions and kaons from D'^ — > K^'k~ 
decays which are selected with high purity without 
using PID information. These calibration events 
are rcweighted in bins of momentum to match 
the momentum distribution of the large J/i/^K^ 
and 'ip{2S)K~^ samples. By this technique, the 
following efficiencies are deduced for DLL/f^ < 6: 
ej^^,^ = (95.8 ± 1.0)%; e^,(25). = (96.6 ± 1.0)%. 
The errors, estimated from simulation, account for 
imperfections in the reweighting and the difference 
of the signal A'+ and 7r+ momenta. 

Partially reconstructed decays populate the 
region below the i?+ mass. _B+/" — > ^K^n 
decays, where the pion is missed, are modelled 
in the kaon-like sample by a fiat PDF with a 
Gaussian edge. A small — > ipK^ir^ component 
is needed to achieve a stable fit. It is modelled 
with the same shape as the partially reconstructed 
5+/° decays except shifted in mass by the B^-B^ 
mass difference, +87 MeV/c^ In the pion-likc 
sample, ■07r"^7r backgrounds are assumed to enter 
with the same PDF, and same proportion relative 
to the signal, as the ipK'^TT background in the 
kaon-like sample. A component of misidcntified 
^ J/i'K+ TT is also included with a fixed 
shape estimated from the data. Lastly, a linear 
polynomial with a negative gradient is used to 
approximate the combinatorial background. The 
slope of this component of the pion-like and 
kaon-like backgrounds can differ. 

The stability of the fit is tested with a large 
sample of pseudo-experiments. Pull distributions 
from these tests are consistent with being normally 
distributed, demonstrating that the fit is stable 
under statistical variations. The yields obtained 
from the signal extraction fit arc shown in Table [2j 

The observables, defined in Eqs. [1] and [2] are 
calculated by the fit, then modified by a set of cor- 
rections taken from simulation. The acceptances 
of ipTT^ and ipK~^ events in the detector are com- 
puted using Pythia [TT] to generate the primary 



Table 2: Raw fitted yields. The labels 'D' and 'U' 
refer to the two polarities of the LHCb dipole. 





B- B+ 


D 

^ U 

u 


528 ± 27 518 ± 27 
421 ± 23 428 ± 23 
13 363 ± 180 13 466 ±181 
10 666 ± 148 11 120 ±155 


D 

K 

U 


94 ± 16 93 ± 16 
82 ± 15 70 ± 13 
2 331 ± 88 2 463± 93 
2 026± 78 1836± 71 



collision and EvtGen [T2] to model the 5+ decay. 
The efficiency of reconstructing and selecting tpir'^ 
and xpK'^ decays is estimated with a bespoke 
simulation of LHCb based on Geant4 pS]- It 
models the interaction of muons and the two 
hadron species with the detector material. The 
total correction e'^^/e'^'^ is 0.985 ± 0.012 and 
1.007 ± 0.021 for i?''/''' and ^"^(25) respectively. 

CP asymmetries are extracted from the ob- 
served charge asymmetries (^Raw) by taking 
account of instrumentation effects. The inter- 
action asymmetry of kaons. A^^^^ is expected 
to be non-zero, especially for low-momentum 
particles. This asymmetry, measured at LHCb 
using a sample of Z?*+ D°tt+ , D° K+n~ 
decays, is —0.010 ± 0.002 if the pion asymmetry 
is zero [14j . The null-asymmetry assumption for 
pions has been verified at LHCb to an accuracy of 
0.25% Hn]. These results are used with enlarged 
uncertainties (0.004, for both kaons and pions) 
to account for the different momentum spectra of 
this sample and those used in the previous analyses. 

In summary, the CP asymmetry is defined as 

^^'=4aw-^Prod-^Lt, (4) 

where the production asymmetry, Apiod, de- 
scribes the different rates with which B~ and B~^ 
hadronise out of the pp collisions. The observed, 
raw charge asymmetry in — s- J/tpK^ is 
-0.012 ± 0.004. Using Eq. g] with the established 
CP asymmetry, A-'/'i'^ = 0.001 ± 0.007 [3], A^.^d 
is estimated to be —0.003 ± 0.009. This is applied 
as a correction to the other modes reported here. 

The different contributions to the systematic 
uncertainties are summarised in Table [3] They 
are assessed by modifying the final selection. 
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Table 3: Summary of systematic uncertainties. The statistical fit errors are included for comparison. 



Simulation uncertainty 


0.045 




0.088 






PID efficiencies 


0.043 




0.052 






j^JI^K (pDG [3]) 




0.0070 




0.0070 


0.0070 


^RaJ^ statistical error 




0.0046 




0.0046 


0.0046 


Detection asymmetries 




0.0056 




0.0056 




Relative trigger efficiency 


0.020 


0.0031 


0.050 


0.0036 


0.0003 


Fixed fit parameters 


0.005 


0.0006 


0.017 


0.0013 


0.0001 


Sum in quadrature (syst.) 


0.065 


0.0106 


0.115 


0.0108 


0.0084 


Fit error (stat.) 


0.110 


0.0268 


0.404 


0.0901 


0.0136 



or altering fixed parameters and rerunning the 
signal yield fit. The maximum variation of each 
observable is taken as their systematic uncertainty. 

The largest uncertainty is due to the use of 
simulation to estimate the acceptance and selec- 
tion efficiencies. It accounts for any bias due to 
imperfect modelling of the detector and its relative 
response to pions and kaons. Another important 
contribution arises from the loose trigger criteria 
that are employed. This uncertainty is estimated 
from the shift in the central values after rerunning 
the fit using only those events where the muons 
passed the software trigger. The use of the PID 
calibration to estimate the efficiency for pions 
to the DLL/fTT < 6 selection also contributes a 
significant systematic uncertainty. 

The measurements of A^'^ depend on the 
estimation of ^prod from the 5+ — > J/ipK^ 
channel. The uncertainty on Apiod is determined 
by the statistical error of A'^^J^ in the fit, the 
uncertainty on the world average of A'^^^'-^ and the 
estimation of A^^^. . These effects are kept separate 
in the table where it is seen that the uncertainty 
on the nominal value of A'^^^-^ dominates. Finally, 
it is noted that the detector asymmetries cancel 
for ^'/'(2S)A' g^j^jj lower systematic uncertainty 
can be reported. 

The measured ratios of branching fractions are 

RJ/'I' = (3.83 ± 0.11 ± 0.07) X IQ-^ 
^V(2S) ^ (3.95 ± 0.40 ± 0.12) X 10"^ 

where the first uncertainty is statistical and the sec- 
ond systematic. R^'^'^^^ is compatible with the one 
existing measurement, (3.99±0.36±0.17) x 10"^ gj_ 
The measurement of R'^/'^ is 3.2(7 lower than the 



current world average, (5.2 ± 0.4) x 10~^ [3]. Us- 
ing the established measurements of the Cabibbo- 
favoured branching fractions [5], we deduce 

B{B^ J/ij-K^) ^ (3.88 ± 0.11 ± 0.15) x 10'^ 
B{B^ i'{2S)T:^) = (2.52 ± 0.26 ± 0.15) x 10~^ 

where the systematic uncertainties are summed in 
quadrature. The measured CP asymmetries, 

A'^P = 0.005 ±0.027 ±0.011 
A^'^^''' = 0.048 ±0.090 ±0.011 
App^^^ = 0.024 ± 0.014 ± 0.008, 

have comparable or better precision than previous 
results, and no evidence of direct CP violation is 
seen. 
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